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Abstract We have studied the effects of thermal 
conduction on the structure of viscous and resistive 
advcction-dominated accretion flows (ADAFs). The 
importance of thermal conduction on hot accretion flow 
is confirmed by observations of hot gas that surrounds 
Sgr A* and a few other nearby galactic nuclei. In this 
research, thermal conduction is studied by a saturated 
form of it, as is appropriated for weakly-coUisional sys- 
tems. It is assumed the viscosity and the magnetic 
diffusivity arc due to turbulence and dissipation in the 
flow. The viscosity also is due to angular momentum 
transport. Here, the magnetic diffusivity and the kine- 
matic viscosity are not constant and vary by position 
and a-prescription is used for them. The govern equa- 
tions on system have been solved by the steady self- 
similar method. The solutions show the radial velocity 
is highly subsonic and the rotational velocity behaves 
sub-Keplerian. The rotational velocity for a specific 
value of the thermal conduction coefficient becomes 
zero. This amount of conductivity strongly depends 
on magnetic pressure fraction, magnetic Prandtl num- 
ber, and viscosity parameter. Comparison of energy 
transport by thermal conduction with the other energy 
mechanisms implies that thermal conduction can be a 
significant energy mechanism in resistive and magne- 
tized ADAFs. This property is confirmed by non-ideal 
magnetoliydrodynamics (MHD) simulations. 
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1 Introduction 

The observational features in active galactic nuclei 
(AGN) and X-ray binaries can be successfully explained 
by the standard geometrically thin, optically thick ac- 
cretion disc model (Shakura & Sunyaev 1973). The mo- 
tion of the matter in the standard thin model of the ac- 
cretion disc is approximately Keplerian, and the energy 
released in the accreting gas is radiated away locally. In 
the past two decades, another type of accretion flow has 
been considered that the energy released due to heating 
processes in the flow may be trapped within accreting 
gas. As, only the small fraction of the energy released 
in the accretion flow is radiated away due to inefficient 
cooling, and most of the energy is stored in the accre- 
tion flow and advected to the central object. This type 
of accretion flow is called as advection-dominated ac- 
cretion flow (ADAF). The models of ADAF have been 
studied by a number of researchers (e.g. Ichimaru 1977; 
Rees et al. 1982; Narayan & Yi 1994; Abramowicz et 
al. 1995; Blandford & Begelman 1999; Ogilvie 1999). 

The observations of black holes confirm existence of 
hot accretion flow that contrasted with classical cold 
and thin accretion disc model (Shakura & Sunyaev 
1973). Hot accretion flows can be seen in supermas- 
sive black holes of galactic nuclei and during quiescent 
of accretion onto stellar-mass black holes in X-ray tran- 
sients (e.g., Lasota et al. 1996; Esin et al. 1997, 2001; 
Narayan et al. 1998a; Menou et al. 1999; Di Mat- 
teo et al. 2000; see Narayan et al. 1998b; Melia & 
Falckc 2001; Narayan 2002; Narayan & Quataert 2005 
for reviews) . Chandra observations provide constraints 
on the density and temperature of gas at or near the 
the Bondi capture radius in Sgr A* and several nearby 
galactic nuclei (Loewenstein et al. 2001; Baganoff et al. 
2003; Di Matteo et al. 2003; Ho et al. 2003). Tanaka & 
Menou (2006) exploited these constraints to calculate 
mean free path of the observed gas. They suggested 
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accretion in such systems are under weakly coUisional 
condition. Moreover, they suggested thermal conduc- 
tion as a possible mechanism by which the sufficient 
extra heating is provided in hot accretion flows. In fol- 
lowing, Johnson & Quataert (2007) studied the effects 
of electron thermal conduction on the properties of a 
spherical hot accretion flows. Their model is applica- 
ble for Sgr A* in the Galactic centre. Because, electron 
heat conduction is important for low accretion rate sys- 
tems. They also found a supervirial temperature in 
the presence of thermal conduction. Similar to Tanaka 
& Menou (2006), they assumed a steady state model, 
but they solved their equations numerically. Abbassi et 
al. (2008) presented a set of self-similar solutions for 
ADAFs with a toroidal magnetic field in which the sat- 
urated thermal conduction has a important role in the 
energy transport in the radial direction. 

Since the observational evidences and magnetohy- 
drodynamics (MHD) simulations have expressed the 
toroidal magnetic field and the magnetic diffusivity 
are important in accretion fiows (see Faghei 2011a and 
references therein), Faghei (2011a) examined the self- 
similar solutions of viscous and resistive ADAFs in the 
presence of a toroidal magnetic field. But, he did not 
consider the effects of thermal conduction in his model. 
While, the recent studies of resistive accretion flows 
have represented that thermal conduction can play an 
importance role in such systems (e.g. Sharma et al. 
2008; Ghanbari et al. 2009). Sharma et al. (2008) 
studied the effects of thermal conduction on magne- 
tized spherical accretion fiows using global axisymmet- 
ric MHD simulations. In their model, when the mag- 
netic energy density becomes comparable to the grav- 
itational potential energy density, the plasma due to 
resistivity is heated to roughly the virial temperature, 
the mean inflow becomes highly subsonic, and most of 
the energy released by accretion is transported to large 
radii by thermal conduction and the accretion rate be- 
came much smaller than Bondi accretion rate. More- 
over, they found that for a larger values of conductive 
parameter, energy transport through thermal conduc- 
tion becomes the dominant energy transport mecha- 
nism at small radii. Ghanbari et al. (2009) studied 
a two-dimensional advective accretion disc bathed in 
the poloidal magnetic field of a central accrctor in the 
presence of thermal conduction. They did not consider 
toroidal component of magnetic field and for simplicity 
assumed the resistivity as a constant. They studied in- 
duction equation of magnetic field in a steady state that 
is not according to anti-dynamo theorem (e.g. Cowling 
1981) and is useful only in particular systems where 
the magnetic dissipation time is much longer than the 
age of the system. On the other hand, this assumption 



implies that the flow is in balance between escape and 
creation of the magnetic field. 

In this paper, we adopt the presented solutions by 
Tanaka & Mcnou (2006) and Faghei (2011a). Thus, 
we will investigate the influences of thermal conduc- 
tion on a viscous and resistive ADAF in the presence 
of a toroidal magnetic field. Moreover, it is assumed 
that magnetic diffusivity in the present model is not 
constant, and escaping and creating of magnetic field 
are unbalanced. From some aspects will be shown that 
the present model is in according with the observations 
and the resistive MHD simulations. The paper is or- 
ganized as follows. In section 2, the basic equations of 
constructing a model for quasi-spherical magnetized ad- 
vection dominated accretion flow with thermal conduc- 
tion will be defined. In section 3, self- similar method 
for solving equations which govern the behaviour of the 
accreting gas was utilized. The summary of the model 
will appear in section 4. 

2 Basic Equations 

Wc suppose a rotating and accreting gas around a 
Schwarzschild black hole of mass M. The fiow is as- 
sumed to be in advection dominated stage, where vis- 
cous and resistive heating are balanced by the advection 
cooling and thermal conduction. We use a spherical 
coordinate (r, 6, </>) centred on the accreting object. 
Furthermore, the fiow is assumed to be steady and ax- 
isymmetric (dt = = 0), and the equations will be 
considered in the equatorial plane, 6 = 7r/2. Thus, all 
flow variables are a function of only r. For the sake 
of simplicity, the general relativistic effects are ignored 
and Newtonian gravity is used. The magnetic fleld in 
the present model has only a toroidal component. Un- 
der assumptions, the model is described by the follow- 
ing equations: 

The continuity equation with mass loss is 

1 d 2 \ 

^"'O = P> (1) 

where p is density, Vr is the radial infall velocity, and p 
the mass-loss rate per unit volume. 
The radial equation of momentum is 

.,^.r,0»-ai-,-iA„a4-li^(.oi). (2, 

where fl is the angular velocity of the fiow, ^Ik = 
^ GM/r^ is the Keplerian angular velocity, Cg the 
isothermal sound speed, which is defined as = Pgas/ Pi 
Pgas being the gas pressure, and ca is Alfven speed, 

which is defined as c\ = B'^I^^^P = ^Vraaglp, Pmag 
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being the magnetic pressure. The angular momentum 
transfer equation is 



LA. 

j,2 ^j. 



vpr 



dr 



(3) 



where the right-hand side of above equation describes 
the effects of viscous torques due to shear {v, here, is 
kinematic coefficient of viscosity). As noted in the in- 
troduction, we assume both of the kinematic coefficient 
of viscosity and the magnetic diffusivity due to turbu- 
lence in the accretion flow. Thus, it is reasonable to 
use these parameters in analogy to the a-prescription 
of Shakura & Sunyaev (1973) for the turbulent. 



11k 



(4) 



where Pm is the magnetic Prandtl number, which is 
assumed a constant of order of unity, rj is the magnetic 
diffusivity, and a is a free parameter less than unity. 
The energy equation becomes 



7 — 1 dr 
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Qdi 



pel d , 2 

(r Vr 



1 r 



2 dr 

Qrad 



(5) 



where Qdiss = Qvis + Qresis is the dissipation rate by 
viscosity Qms and resistivity Qresis, Qrad represents 
the energy loss through radiative cooling, and Qcond is 
the energy transported by thermal conduction. For the 
right-hand side of the energy equation, we can write 



Cadv 



diss 



Qrad ^" Q 



cond 



where Qadv is the advective transport of energy. We 
exploit the advection factor, / = 1 — Qrad/ Qdiss, that 
describes the fraction of the dissipation energy which 
is stored in the accretion flow and advccted into the 
central object rather than being radiated away. In gen- 
eral, the advection factor depends on the details of the 
heating and radiative cooling mechanisms and will vary 
by position (e.g. Watari 2006, 2007; Sinha et al. 2009). 
Here, we assume a constant / for simplicity. Clearly, 
the flow in the case of / = 1 is in the extreme limit of no 
radiative cooling and in the limit of efficient radiative 
cooling, we have / = 0. 

As mentioned, the inner regions of hot accretion 
flows are collisionless and the electron mean free path 
due to Coulomb collision is larger than the radius. This 
property is described as saturation (Cowie & McKee 
1977). The traditional equation of heat flux due to 
thermal conduction, Fcond = —k\JT which k being the 
thermal conduction coefficient, is not valid for such sys- 
tems. Because this equation is suitable for the colli- 
sional plasma, in which mean free path of electron en- 
ergy exchange is smaller than temperature scale height. 



Cowie & McKee (1977) derived the heat flux for the col- 
lisionless plasma as 



(7) 



where 0^ is a factor is less than unity and is called 
as saturation constant. Now, the viscous and resistive 
heating rates and the energy transport by thermal con- 
duction are expressed as 



vpr" 



dr 



47r 



J.2 Qj, 



(8) 



(9) 



(10) 



where J = V x B is the current density, B being the 
magnetic field. We used a minus sign in equation (10). 
Because we want thermal conduction as energy trans- 
port mechanism outward. On the other hand, heat con- 
duction flux will be behaved like a cooling mechanism 
in accretion flow. Finally, since we consider only the 
toroidal field, the induction equation with field escape 
can be written as 



r dr 



rVrBu, - ri—{rB^) 
dr 



(11) 



(6) where B^ is the toroidal component of magnetic field 



and B^p is the field escaping/creating rate due to a mag- 
netic instability or dynamo effect. This induction equa- 
tion is rewritten as 



ld_ 

r dr 



A-Kpc 



a 



1 



4/3P™ rpVlK dr 



(12) 



where /3 is the degree of magnetic pressure to the gas 
pressure and can be defined by 



1 



Pgas 



(13) 



In this paper, we apply the steady self-similar methods 
to solve the system equations. Thus, this parameter 
will be constant throughout the disc. While, Khesali 
& Faghei (2008, 2009) showed that it varies by posi- 
tion. In hot accretion flows, typical value of /3 is in the 
range 0.01-1 (e.g. De Vflliers et al. 2003; Beckwith 
et al. 2008). Here, we will also consider the magneti- 
cally dominated case (/3 > 1). Because, when thermal 
instability happens in an ADAF, the MHD numerical 
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simulations imply that the thermal pressure rapidly de- equations [(l)-(3), (5), and (11)], we can obtain the 



creases while the magnetic pressure increases due to the 
conservation of magnetic flux (Machida et al. 2006). 
This will result in large /3 and forms a magnetically 
dominated accretion flow (Bu et al. 2009). 



following relations: 



s + 



apoci 



(21) 



3 Self-Similar Solutions 



3.1 Analysis 



C2 - 1 - C3 [S - 1 + /?(1 + S)] , 



(22) 



The self-similar method is useful to understand of the 
physics of accretion flows. This method is familiar due 
to its wide range of applications in many research fields 
of astrophysics. Self-similar solution, although consti- 
tuting only a limited part of problem, is often useful to 
understand the basic behaviour of the system. Thus, 
in order to seek similarity solutions for the above equa- 
tions, we seek solutions in the following form: 



ci = 3(s + 2)c3, 



aci 



2{s - 1) + 37 



7-1 



(23) 



-lOfZ-sCg/' 



1 

^+2 



(24) 



GA'L 



(14) 



(15) 



Bo = -—GM,^27rpo/3c3 



2ci + |L(i + ,) 



(25) 



Above equations express for s = —3/2, there is no mass 
loss, while for s > —3/2 mass loss (wind) exists. After 
algebraic manipulations, we obtain an algebraic equa- 
tion for C3: 



GAL 



Anp r 



(16) 



(17) 



where coefficients ci, C2, and C3 are determined later. 
We assume a power-law relation for density 



p{r) = por' 



(18) 



where pq and s arc constant. By using above self- 
similar quantities, the mass-loss rate and the field es- 
caping/creating rate must have the following form: 



9af ( 2/3 



3a(s-|-2) 2(s-l) + 37 
2 ^ 7- 1 



9P„ 



{s + iy + {s + l)l3 + s-l 



V^-- /a = 0,(26) 



and the rest of the physical variables are 



po = -SapoVGAhis + 2) (^s + 



C3, 



ci = 3c3(s -t- 2), 



(27) 



(28) 



p{r) = por' 



-3/2 



(19) 



l-^(,s + 2)2c2+c3 [(,s + l)/3 + s 



(29) 



(20) 



where po and Bq are constant. 

Substituting the above solutions in the continuity, 
momentum, angular momentum, energy, and induction 



-3asGAL{s + 2)^/^03/ VSttpo^ x 
s + l 



6P™(s + 2) 



(30) 



We can solve algebraic equation (26) numerically and 
clearly only real roots which correspond to positive ci 
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Fig. 1 Physical quantities of the flow as a function of Fig. 2 Same as Fig. 1, but /? = 1.0, and solid, dotted, 

saturation constant for 7 = 1.3, a — 0.5, / = 1, and dashed, and dot-dashed lines represent Pm = cxd, 1, 2/3, and 

Pm = 100. Solid, dotted, dashed, and dot-dashed lines rep- 1/2. 
resent /3 = 0.1, 0.3, 0.7, and 1.5. 



are physically acceptable. Without thermal conduc- 
tion, i.e. (ps = 0, (26) and similarity solutions reduce 
to Faghei (2011a). But our main algebraic equation 
includes thermal conduction. 

3.2 Numerical Results 

Now, we consider the behaviour of solutions in the pres- 
ence of thermal conduction. But in this paper, only case 
of no wind (s = —3/2) is considered that po = and 
B^p cx r~^^/^. In addition to introduced coefficients, we 
also define a new parameter of C4 that is the right-hand 
side of equation (24) 



C4 = af 
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+ -^C3(l + s)^ 



10(j)sC- 



1/2 



(31) 



Equations (5), (24), and (31) imply that the parameter 
C4 is the advection transport of energy. The behaviour 
of coefficients Ci as a function of (j)s are shown in Figures 
1-3. Moreover, Figure 1 represents the profiles of physi- 
cal quantities for several values of the magnetic pressure 
fraction, i. e. f3 = 0.1,0.3,0.7, and 1.5. The value of /3 
measures the strength of magnetic field, and a larger (3 
denotes a stronger magnetic field. Figure 2 represents 
the profiles of physical quantities for several values of 



magnetic Prandtl number, i. e. Pm = 00, 1,2/3, and 
1/2. The smaller values of Pm denotes a stronger mag- 
netic diffusivity, 77. Figure 3 shows the profiles of phys- 
ical quantities for several values of adiabatic index, i. 
e. 7= 1.2,1.25,1.3, and 1.35. 

The solutions in Figures 1-3 imply that the radial 
infall velocity, ci, and the sound speed, C3, both de- 
crease with the magnitude of conduction, while the 
squared angular velocity, c|, increases. These proper- 
ties are qualitatively consistent with dynamical analysis 
of Faghei (2011b). One and two dimensional simula- 
tions of hot accretion flows have also shown that ther- 
mal conduction reduces the flow temperature (Sharmal 
et al. 2008; Wu et al. 2010). The profiles of advec- 
tion transport of energy, C4, in Figures 1-3 imply that 
thermal conduction behaves as a cooling mechanism, 
resulting in a local decrease of the gas temperature rel- 
ative to the original ADAF solution. At the same time, 
the gas adjust its angular velocity (which increases the 
level of viscous dissipation) and reduces its infiow speed. 

Figure 1 shows the radial velocity, sound speed, and 
advection transport of energy decrease by adding the 
magnetic pressure fraction, /3. These properties are 
qualitatively consistent with results of Bu et al. (2009) 
and Faghei (2011a). Moreover, the angular velocity de- 
creases with the magnitude of magnetic field that is in 
according with results of Khesali & Faghei (2009) and 
Faghei (2011a). 
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Fig. 3 Same as Fig. 1, but /3 = 1.0, and solid, dotted, 
dasiied, and dot-daslied lines represent 7 = 1.2, 1.25, 1.3, 
and 1.35. 



Figure 2 shows the magnetic diffusivity has the oppo- 
site effects of thermal conduction on the physical vari- 
ables. As, by adding the magnetic diffusivity, P,^^, the 
radial velocity, sound speed, and advection transport of 
energy increase, while the rotational velocity decrease. 
These results are similar to resistive ADAF models 
without thermal conduction (e. g. Faghei 2011a). 

Figure 3 represents the gas adiabatic index similar 
to magnetic diffusivity has the opposite effects of ther- 
mal conduction on the physical quantities. As, with 
the magnitude of 7, the inflow and sound speed in- 
crease, while the rotational velocity decreases. This is 
in accord with dynamical study of hot accretion flow 
(Faghei 2011b). Moreover, Figure 3 shows that the gas 
adiabatic index contributes with thermal conduction to 
reduce advection transport of energy. It can be due 
to decrease of rotational velocity by adding 7, which 
reduces the level of viscous dissipation. 

The studies of hot accretion flows (e. g. Shadmchri 
2008) imply that the solution for a given set of the input 
parameters reaches to a non-rotating limit at a specific 
of 0s which we denote it by 0^. With zero insertion of 
C2 in equations (22)-(24), 0^ can be written as 



1 

60 



/3/ 5/3-7 ' 

7-1 , 



2(/3 + 5) 




(32) 



We can not extend the studies beyond 0^, because the 
right-hand side of equation (29) becomes negative and 
a negative is clearly unphysical. As, the rotational 
velocity profiles in Figures 1-3 show, we have selected 
the input parameters that is positive. The behaviour 
of critical saturation constant, 0^, as a function of /3 for 
several values of magnetic Prandtl number is shown in 
Figures 4. In left panel of Figure 4, the viscosity param- 
eter value is a = 0.1, and in right panel is a. = 0.2. The 
profiles of Figure 4 show the critical saturation constant 
highly depends on magnetic pressure fraction, /3, mag- 
netic Prandtl number, P^, and viscosity parameter, a. 
As, higher values of /3 corresponds to larger 0^. Critical 
saturation constant also increases with higher value of 
a. Since, magnetic diffusivity is proportional to P,^^, 
Figure 4 shows that the magnetic diffusivity similar to 
magnetic field increases (jf^ value. 

As mentioned in the introduction, Sharma et al. 
(2008) by resistive MHD simulation studied a spher- 
ical accretion with thermal conduction. They found 
for even modest thermal conductivities, conduction is 
the significant mechanism of energy. Here, to compare 
conduction mechanism to others, we study the ratio 
of energy transport by thermal conduction, Qcond, to 
the gas heating rate by viscosity Qvis and resistivity 
Qresis- Such solutious are shown in Figure 5 for two 
cases of non-resistive (left-panel) and resistive (right- 
panel) flows. The solutions imply that thermal con- 
duction is the significant energy mechanism in the flow. 
This result confirms simulation of Sharma ct al. (2008). 
Moreover, the ratio of Qcond to Qdiss increases slightly 
by adding the magnetic field and docs not change for 
different values of magnetic Prandtl number. Because, 
a large fraction of Qdiss is generated by viscous dissi- 
pation. 



4 Summary and Discussion 

The collision timescale between ions and electrons in 
hot accretion flows is longer than the inflow timescale. 
Thus, the inflow plasma is coUisionless, and transfer 
of energy by thermal conduction can be dynamically 
important. The low coUisional rate of the gas is con- 
firmed by direct observation, particularly in the case of 
the Galactic centre (Quatacrt 2004; Tanaka & Mcnou 
2006) and in the intraclustcr medium of galaxy clusters 
(Sarazin 1986). 

In this paper, the structure of a magnetized ADAF in 
the presence of resistivity and thermal conduction is in- 
vestigated. We assumed the magnetic field has a purely 
toroidal component. We adopted the presented solu- 
tions by Tanaka & Menou (2006) and Faghei (2011a). 
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a = 0.1 a = 0.2 P^'" -P„ = 0-5 




Fig. 4 The critical saturation constant as a function of 
the ratio of magnetic pressure to gas pressure. Solid, dotted, 
dashed, and dot-dashed lines represent Pm = 10, 5, 1, and 
0.5. The input parameters are set to 7 = 5/3, / = 1, 
s = —3/2, and the viscous parameter a in left-panel is 0.1 
and in right-panel is 0.2. 



Fig. 5 The ratio of energy transport by thermal con- 
duction, Qcond, to the gas heating rate by viscosity Qvis 
and resistivity Qresis as a function of saturation constant. 
Solid, dotted, dashed, and dot-dashed lines represent /3 — 
0.0,1.0,3.0, and 5.0. The input parameters are set to 
7 = 4/3, f = 1, s — —3/2, a = 0.5, and the magnetic 
Prandtl number in left-panel is 00 and in right-panel is 0.5. 



Thus, we assumed that angular momentum transport 
is due to viscous turbulence and the a- prescription is 
used for the kinematic coefficient of viscosity. We also 
assumed the flow does not have a good cooling effi- 
ciency and so a fraction of energy accretes along with 
matter on to the central object. In order to solve the 
equations that govern the structure behaviour of mag- 
netized ADAF with thermal conduction, we have used 
steady self-similar solution. 

The solutions showed the radial infall velocity and 
sound speed in the presence of thermal conduction both 
decrease, while angular velocity increase. These prop- 
erties are consistent with dynamical study of hot ac- 
cretion flow (e. g. Faghei 2011b) and from some as- 
pects also are in accord with simulations of Sharma et 
al. (2008) and Wu et al. (2010). Moreover, the so- 
lutions represent the magnetic diffusivity and thermal 
conduction have the opposite effects on physical quanti- 
ties. For a moderate thermal conduction, the solutions 
imply that thermal conduction can play an important 
role in energy mechanism of the system. This property 
is qualitatively consistent with non-ideal simulations of 
Sharma et al. (2008). 

In the present model, accretion flow is studied in 
one-dimensional approach and ignored from latitudinal 
dependence of physical quantities. There are some re- 
searches in two-dimensional approach that express the 
importance of such studies (Tanaka & Menou 2006; 
Ghanbari et al. 2009; Wu et al. 2010 ). Thus, latitudi- 
nal study of present model can be investigated in other 
research. Here, we used a saturated heat conduction 
flux. While, there are some studies with unsaturated 
heat flux (e. g. Shcherbakkov & baganoff 2010) that 
show a good agreement with observations. Thus, the 



study of the present model in a unsaturated case will 
be interesting. 
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